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The new definition of lattice gauge fields and the Landau gauge 
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The Landau gauge fixing algorithm in the new definition of gauge fields is presented. In this algorithm a new 
solver of the Poisson equations based on the Green's function method is used. Its numerical performance of the 
gauge fixing algorithm is presented. Performance of the smeared gauge fixing in SU(3) is also investigated. 



1. The Landau gauge fixing algorithm in 
the new definition of gauge fields 



We define the gauge field on links as an 
element of SU (3) Lie algebra as, 



U 



where Al 



-A 



1) 



Our definition is different from usual conven- 
tion[p| in [/-linear form, 



AW = ±(U-U%, 



acclcss part 



(2) 



and it is more natural as an averaged gauge field 
on the link in the continuum theory. We use the 
analytic method of eigenspace projection opera- 
tor together with Cardano's eigenvalue solver for 
transformation between U and A. It is remark- 
able that the definition (|l|) of gauge field allows 
the same type of formulation of the lattice Landau 
or Coulomb gauge condition as the continuum 
theory, i.e., minimizing condition of \\A\\ 2 along 
the gauge orbits H(§. The definion A — \ogU re- 
quires a proper choice of branch for traceless A, 
and A becomes a nonanalytic function of U, and 
yet || A || 2 still remains a continuous function of U. 
This fact justifies the existence of the global mini- 
mum of || A 5 1| 2 along the gauge orbit (the minimal 
Landau gauge condition). 
The gauge transformation is 



A3 +4 
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(3) 



where g x = e £x and e x is a traceless antihermitian 
matrix. 

When e is infinitesimal 

SA x>fl = {T>e) x>fl = f(adj Ax ,Jde x>fl + [A x>ll , e x+ll/2 ] 
x/2 



where 



and 



tanh(a;/2) 
_ 1 f \ 



(4) 



Gribov region ft HI is defined as a set of local 
minimum points on the gauge orbit F 

A||/ll 2 = -2{6M|e) + (e| - 3D\e) + ■ ■■ (5) 

ft = {A\ - dV > , dA = 0} . (6) 

The minimal Landau gauge condition de- 
fines the fundamental modular region by the 
abosolutc minimum along the gauge orbit, 



A = {A\\\A\\ 2 = Min s ||A 9 || 2 }, Ac ft 



(7) 



The equation to solve for the Landau gauge 
fixing is highly nonlinear, 



dA 9 — where g = e e 



(8) 



Linearizing this equation, we can solve it iter- 
atively by the Newton method with an over- 
relaxation factor 77, 



e = (-dViA))- 1 !! dA 



(9) 



This is a non-Abelian extension of the Fourier 
acceleration type|]]||. We make a perturbation 
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expansion with respect to A for the inverse of the 
covariant Laplacian, as 

i-dviA))- 1 e ((-d'r'dDn-d 2 )- 1 (io) 

n=0 

where V = T>(A) - d. 

Performance check for various N enc t is given be- 
low. 

2. The new solver of the Poisson equation 

The inversion of the Laplacian is an important 
problem in numerical analyses, and performances 
of various methods have been investigated^. Ef- 
ficiency of the methods depends on types of the 
problem posed, i.e., balance of accuracy and CPU 
time. We report a new trial, use of the Green 
function. While high accuracy being maintained, 
direct use of the Green function in a single ma- 
chine computaion requires a work of order V 2 , 
and is obviously inefficient in larger size V. We 
devise a method that makes a partial use of the 
Green function for coarse lattice with a precondi- 
tioning for the Poisson equation. 

Let us consider d dimesional Poisson equation, 

-d 2 4> = 2d(I-A)4> = p, (11) 

where A stands for an averaging operation on 
neighboring sites. A preconditioned equaiton is 
given as 

2d(I-A 2 )<p= (I + A)p, (12) 

where (/ — A 2 ) becomes block-diagonal on 
even(odd)-site. After solving the even-site equa- 
tion (p"2|), we use ( p"l| ) for 4> on odd-site. 

The Green function for even(odd)-site equation 
( p"l"| ) is defined as 

2d(I A 2 )G = 8 . (is) 

we found that a good approximant G of the Green 
function G can be given by the Green function G c 
on coarser lattice with volume V/2 d , as 

G = (a + a 1 A 2 +a 2 A A )5-c + bA (i G c , (14) 

where ao, ffli, a,2, b and c are optimized param- 
eters. We can add an extra function AG for 



a higher accuracy correction. For an accuracy, 
10 , of the Green function G on 8 3 x 16 lattice, 
numbers of sites for supports of G c and AG are 
256 and 162, respectively, out of 4096 of full even- 
site. Owing to the high accuracy of the above 
Green function algorithm, this Poisson equation 
solver allows the over-relaxation factor rj to be a 
larger value, 1.5 to 1.8 in gauge fixing algorithm 
(||) and (|l^) and gives a 30 % faster performance 
than that in the usual conjugate gradient solver. 

The following Fig. 1 shows processs of gauge 
fixing in use of the Green function method and 
the conjugate gradient method. 
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Figure 1. A comparison of ZogioMax|cL4| in 
the conjugate gradient method (dashed) and the 
Green's function method(continuous). 



3. Numerical performance of the gauge fix- 
ing algorithm 

We compare performance of the gauge fixing 
algorithm @ and @ for various A perturbation 
orders N en d in Fig. 2, which shows that the non- 
Abelian extension improves efficiency. Relative 
CPU time ratios over that of N en d = 3 case are 
1.05 for N en d — 4 case and 1.91 for N en d = 2 
case. 
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Figure 2. The dependence of logioMax\dA\ on 
the inclusion of the 2nd(dashed), 3rd(continuous) 
and 4th(dash-dotted) order perturbative commu- 
tator term in the Green's function method. 



Figure 3. A histogram of the number of Gribov 
copies obtained after the random gauge transfor- 
mation. Transformation of the 50 samples to the 
unique minimum ||^4|| 2 = 0.2430 was achieved. 



In general, the Fourier accelaration type of 
methods accomplish much better performance in 
ultimate precision over simple norm-minimizing 
methods, and our method is free of critical slow- 
ing down. 

4. Smeared gauge fixing in SU(3) 

As is well known, the Landau gauge fixing suf- 
fers from the problem of Gribov's copies, that is, 
the configuration A 9 is captured by the plenty 
of local minima of ||A ff || 2 along the gauge orbits. 
In order to avoid noise of those Gribov's copies 
when measuring, e.g., the gluon propagators etc., 
the method of gauge fixing to the global mini- 
mum of the 1 1 A 9 1| 2 is necessary. Hetrick and de 
Forcrand proposed a clever method of smeared 
gauge fixing^. They investigated performance 
of the smeared gauge fixing in case of SU(2), and 
reported that it works well for large (3. Here we 
report that it works very well also for SU(3) in 
case of large (3. We observed that it works per- 
fectly in gauge fixing to the global minimum of 
|| A 9 1| 2 , by using 50 configurations generated by 
random gauge transformation from a sample in 



(3 = 5 on 4 3 x 8 lattice, all of which are trans- 
formed to the unique minimum. 

REFERENCES 

1. F. Di Renzo, Nucl. Phys. B(Proc. Suppl.)53 

(1997) 819. 

2. H. Nakajima and S. Furui, talk 
in YKIS97(1997); idem, talk in Confinement 
111(1998); idem, in this poceedings. 

3. T. Maskawa, H. Nakajima, Prog. Theor. 
Phys. 60 (1978) 1526, Prog. Theor. Phys. 63 
(1980) 641. 

4. M. A. Semenov-Tyan-Shanskii, V. A. Franke, 
Zapiski Nauchnykh Seminarov, Leningrad- 
skage Otdelleniya Mathemat. Instituta im 
Steklov ANSSSR 120 (1982) 159. 

5. V.N. Gribov, Nucl. Phys. B139 (1978) 1. 

6. D. Zwanziger, Nucl. Phys. B364 (1991) 127. 

7. C.T.H. Davies et al., Phys. Rev. D 37 (1988) 
1581. 

8. A. Cucchieri and T. Mendes Phys. Rev. D 57 

(1998) R3822. 

9. J.E. Hetrick and Ph. de Forcrand, Nucl. Phys. 
B(Proc. Suppl.)63(1998) 838. 



